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ABSTRACT 

Previous work bv iMotch et alJ Jl985h suggested that in the low/ha rd state of GX 339-4 the soft X-ray power- 
law extrapolated backward in energy agrees with the IR flux level. ICorbel & Fended J2002I) later showed that 
the typical hard state radio power-law extrapolated forward in energy meets the backward extrapolated X-ray 
power-law at an IR spectral break, which was explicitly observed twice in GX 339-4. This TR coincidence’ 
has been cited as further evidence that synchrotron radiation from a jet might make a significant contribution 
to the observed X-rays in hard state black hole systems. We quantitatively explore this hypothesis with a series 
of simultaneous radio/X-ray observations of GX 339-4, taken during its 1997, 1999, and 2002 hard states. 
We fit these spectra, in detector space, with a simple, but remarkably successful, doubly broken power-law 
model that indeed requires an IR spectral break. For these observations, the break position and the integrated 
radio/IR flux have stronger dependences upon the X-ray flux than the simplest jet model predictions. If one 
allows for a softening of the X-ray power law with increasing flux, then the jet model can agree with the 
observed correlation. We also And evidence that the radio flux/X-ray flux correlation previously observed in 
the 1997 and 1999 GX 339-4 hard states shows a ‘parallel track’ for the 2002 hard state. The slope of the 
2002 correlation is consistent with observations taken in prior hard states; however, the radio amplitude is 
reduced. We then examine the radio flux/X-ray flux correlation in Cyg X-1 through the use of 15 GHz radio 
data, obtained with the Ryle radio telescope, and Rossi X-ray Timing Explorer data, from the All Sky Monitor 
and pointed observations. We again And evidence of ‘parallel tracks’, and here they are associated with ‘failed 
transitions’ to, or the beginning of a transition to, the soft state. We also And that for Cyg X-1 the radio flux is 
more fundamentally correlated with the hard, rather than the soft. X-ray flux. 

Subject headings: accretion, accretion disks - black hole physics - radiation mechanismsmon-thermal - X- 
rays:binaries 


1. INTRODUCTION 

Both Cyg X-1 and GX 339-4 in their spectrally hard, radio- 
loud states have served as canonical examples of the so-called 
‘low state’ (or ‘hard state’) of galactic black hole candidates 
(see iPottschmidt et alJ 120031 iNowak et alJ l2003. and refer¬ 
ences therein). As we discuss below, in this state the X- 
ray spectrum is reasonably well-approximated by a power- 
law with photon spectral index (i.e., photon flux per unity 
energy oc E~'^) of T w 1.7, with the power-law being ex¬ 
ponentially cutoff at high energies (« 100keV). It long has 
been suggested that such spectra are due to Comptoniza- 
tion of soft photons from an accretion disk by a hot corona 
in the central regions of the compact object system (e.g., 
ISunvaev & Trumped Il979h . Comptonization models have 
been very successful in describing t he broad-band X-ray/sof t 
gamma-ray spectra of both Cyg X-1 jPottschmidt et al.l2003h 
and GX 339-4 JNowak et alJ2002^ . 

It recently has been hypothesized, however, that the X-ray 
spectra of hard state sources might instead be due to syn¬ 
chrotron and synchrotron self-Compton (SSC) radiation from 
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a mildly relativistic jet dMarkoff et al.ll2fl(fill2flfl3h . Jet mod¬ 
els have been prompted in part by multi-wavelength (radio, 
optical, X-ray) observations of hard state systems. Although 
radio emi ssion from Cyg X-1 w as first observed quite some 
time ago JBraes & Milevi ll 971ll. it is more rec ently that a ra¬ 
dio jet has been imaged 7strrling et alJl200i1) . and that the 
low/hard state X-ray flux h as been shown to be positively cor¬ 
related with the radio flux jPoolev et al.l1998l) . 

Radio emission also h a s been discovered in GX 339-4 
llSood & Camnell-Wilsonl Il994li . The radio emission 
is correlated with X-ra y flux in spectrally hard states 
(iHannikainen et al l 1 1998ft. but is q uenched during spec- 
trallv soft states Jpender et al.lll999l) . Furthermore, in hard 
states of GX 339-4, the 3-9 keV X-ray flux (in units of 
10“*** ergcm“^ s“*) is relat ed to the 8.6GHz r adio flux (in 
mJy) by Ex ~ 0.46/E* JCorbel et alJ 120031 noting that 
throughout this paper we shall use caligraphic script to denote 
flux densities, i.e., flux per unit energy, and roman script to 
denote flux integrated over an energy band). This correlation 
was seen to hold over several decades in X-ray flux, and also 
to hold for two hard state epochs that were separated by a pro¬ 
longed, intervening soft state outburst. Similar correlations 
were found bet ween energy bands in the 9-200 ke V range and 
the radio flux JCorbel et alJl2003l) . It further has been sug¬ 
gested that the Ex oc /E/ correlation is a universal property of 
the low/hard state of black hole binaries llGallo et alJ2003ft . 

This specific power-law dependence of the radio flux upon 
the X-ray flux naturally arises in synchrotron jet models 
jpalcke & Riermannll199 5UCorbel et ^l2003t iMarkoff et alJ 

I2003I lHein7&'Sunvaen2003ft . Both the location of the 
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Fig. 1.— Left: Unfolded spectra of a simultaneous radio and X-ray spectrum of GX 339—4 (D, Obsei'vation ID 40031-03-01, in Tables 1 and 2), fit with an 
absorbed, exponentially cutoff, doubly broken power-law and a gaussian line. Residuals are from the proper forwai'd folded model fit. Right: Same unfolded 
spectrum as on the left, showing just the RXTE data. 


break from an optically thick* to an optically thin radio 
spectrum (presumed to continue all the way through the X- 
ray), and the amplitude of the optically thick portion of 
the radio spectrum, vary with input power to the jet so 
as to produce the Fx oc scaling. Assuming, however, 
that the X-ray power of a disk/corona is proportional to 
M^, where M is the accretion rate, while the jet power is 
proportional to M, a l so reproduces the scalin g relationship 

iFender et alJ 

120031 iMerloni et alj2003tlFalcke et al j2004. _ 

Interestingly, nearly 20 years ago lMotch^ al.l Jl985h noted 
that for a set of simultaneous IR, optical, and X-ray observa¬ 
tions of the GX 339-4 hard state, the extrapolation of the X- 
ray power-law to low en ergy agreed with the over all flux level 
of the optical/IR data. ICorbel & Fended J2002h reanalyzed 
these observations, which did not include simultaneous radio 
data, as well as a set of (not strictly simultaneous) radio/IR/X- 
ray observations from the 1997 GX 339-4 hard state. They 
showed that the low energy extrapolation of the X-ray power- 
laws, and the high energy extrapolation of the radio power- 
law, coincided with a spectral break in the IR. The spectral 
shapes below and above the IR break were roughly consis¬ 
tent with the radio and X-ray power-laws, respectively. This 
TR coincidence’ has been cited as further evidence that the 
jet not only accounts for the flat spectrum from radio through 
infrared, but also that optically thin emission from the jet, oc¬ 
curring at energies above the IR brea k, provides a significant 
contribution to the observed X-rays ([Corbel & Fende3l2002t 

ICorbel et alJ200llMarkoff et alJ200T 

In this paper, we quantitatively examine the TR coinci¬ 
dence’ with a series of simultaneous radio/X-ray spectra of 
GX 339-4. We use broken power-law hts to the combined 
radio and X-ray spectra to determine the extrapolated posi¬ 
tion of the break as a function of observed X-ray flux. We 
then reassess the correlation between radio and X-ray flux for 
GX 339-4 by including more recent hard state observations 
that occurred at fairly high X-ray fluxes. Based upon our re¬ 
sults for GX 339-4, we also reassess this correlation for the 
case of Cyg X-1. We then summarize our results. 

* By optically thick, we mean radio energy flux density Si, oc 1 '°’’, with 
Or > 0. Throughout this work, for both radio and X-ray spectra, we will 
follow the convention o = 1 — F. 


2. THE TR COINCIDENCE’ IN GX 339-4 
2.1. Data Analysis 

We consider a set of ten simultaneous radio/X-ray observa¬ 
tions o f GX 339-4, eight of which were discussed previously 
by us JWilms et al.l 119991 iNowak et alJ I2002I ICorbel et alJ 
120031) and come from the 1997 or 1999 hard s tate, and two of 
which were discussed bv iHoman et lil ilTtiol and come from 
the 2002 hard state (approximately a month before a soft state 
transition). All X-ray observations were performed with the 
Rossi X-ray Timing Explorer (RXTE). Their observation IDs 
and associated radio flux densities and integrated X-ray fluxes 
are presented in Table 1. Note that four of these observations 
are further labeled A-D, as we single these out for special 
discussion below. Observations A and B (40108-01-03 and 
40108-01-04) occurr ed immediately afte r the 1999 soft-to- 
hard state transition JNowak et alJl2003) and have optically 
thin radio spectra (or < 0). Observation C and D (70109- 
01-02 and 40031-02-01) have the brightest X-ray fluxes in 
our sample, and are among the brightest hard X-ray states ob¬ 
served in GX 339-4 to date. 

To analyze the X-ray spectra of these observations, RXTE 
response matrices were created using the software tools avail¬ 
able in HEASOFT 5.3^. The Proportional Counter Array 
(PCA ) data were rebinned to have a minimum of 30 counts 
per bin, uniform systematic uncertainties of 0.5% were ap¬ 
plied, and only data between 3 and 22 keV were considered. 
The Eligh Energy X-ray Timing Explorer (EIEXTE) data were 
coadded from the two individual clusters and then were re¬ 
binned to have a minimum signal-to-noise ratio (after back¬ 
ground subtraction) of 10 in each bin. We considered EIEXTE 
data only in the 18-200keV range. In the fits discussed be¬ 
low, a multiplicative constant was allowed between the PCA 
and EIEXTE normalization, with the PCA constant fixed to 
unity. The normalizations of the two instruments were always 

^ The use of HEASOFT 5.3 is very important here, as we find extremely 
good agreement between the Proportional Counter Array and High Energy X- 
ray Timing Explorer when fitting power-law models to the Crab pulsar plus 
nebula system. This is true for both the power-law normalization and slope, 
both of which must be determined very accurately when extrapolating over a 
large range of energies between the radio and X-ray spectra. This spectral and 
flux agreement is in marke d contrast to eai'lier versions of HEASOFT (see, for 
example, the discussion of lWilms et alJ199^ . 
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Fig. 2.— Unfolded RXTE spectra of Cyg X-1 (Wilms et al., in prep.), 
fit with an absorbed, exponentially cutoff, broken power-law and a gaussian 
line. Residuals are from the proper forward folded model fit. 


found to be the same to within a few percent (see Table 2). 

The radio data for observation D (40031-03-01) were ob¬ 
tained with the Australia Telescope Compact Array (ATCA) 
at 4.8 GHz and 8.6 GHz. The radio data for observation C 
(70109-01-02) were also obtained with ATCA, but at 1.4 and 
2.4 GHz. These radio data, as well as the data for observation 
B (40108-01-04), the 1999 April 22 observation (40108-02- 
02), and the 1999 May 14 observation (40108-02-03) have 
been reanalyzed by us. All oth e r radi o data come from 
I Wilms et all (119991) . iNowak et al.l (l2002ft . and ICorbel et al.l 
(120001) . In cases of small discrepancies between these lat- 
ter references, we have adopted the radio flux values of 
ICorbel etal.l l200g). since, aside from the radio observations 
reanalyzed by us, that work represents the most up to date 
analysis of the radio data presented here. 

Note that we have considered only data wherein we have 
nearly simultaneous radio an d X-ray spect r a. Th e radio/X- 
ray observations discussed bv iNowak et aDdlOOl all have a 
large degree of overlap between the radio and X-ray observ¬ 
ing times, and they show very little long time scale (> 1000 s) 
variability (ICorbel et al.l200dh . Radio observations C & D oc¬ 
cur approximately a half day after their corresponding RXTE 
observations, but, again, neither the radio data nor the RXTE 
data show any evidence for long term variability within the 
observations themselves. Still, the lack of strict simultane¬ 
ity might be a relevant factor in some of the results discussed 
below. 

We further require that the source was bright enough to fit 
the X-ray spectru m in both the PCA and t he HEXTE data. 
As dis cussed by ICorbel et alJ (120031) and IWardzihski et alJ 
Iliol . at low flux levels there is background contamination 
that affects the RXTE spectra of GX 339-4. The spectra 
can be corrected adequatel y to yield reasonab ly accurate in¬ 
tegrated soft X-ray fluxes ICorbel et al.ll2003l) : however, the 
contamination is much more problematic for the spectral ex¬ 
trapolations discussed here. 

The observations were analyzed with the Interactive 
Spectral Interpret ation System (ISIS; 

iHouck & Denicol^ l2000l) . For our purposes, there are 
three major reasons for our use of ISIS. First, ISIS 
uses s-lang as a scripting language and hence has most 
of the programmability of IDL (Interactive Data 
Language) or MATLAB, while retaining all the models 


of XSPEC (lArnaudI Il996h . Second, data input without 
a response matrix (i.e., the radio data) are automatically 
presumed to have an associated diagonal response with one 
cm^ effective area and one second integration time. Thus, we 
used a fairly straightforward s-lang script to convert the 
radio data from mJy to photon rate in narrow bands around 
the observation frequencies. This then was used as input 
for the simultaneous radio/X-ray fits. (So long as they are 
narrow, the widths of the radio bands do not affect the fits.) 
We then set the fractional error bars of these ‘count rate’ data 
equal to the fractional error from the radio measurements. 

The third reason for using ISIS is that it treats ‘unfolded 
spectra’ (shown in Fig. Q in a model-independent manner. 
The unfolded spectrum in an energy bin denoted by h, as used 
to create Fig.^ is defined by: 


Funfold (^) — 


lC(h)-B(h)]/At 
f R(h,E)A(E)dE ’ 


( 1 ) 


where C(h) is the total detected counts, B{h) is the background 
counts, Af is the integrated observation time, R{h,E) is the 
unit normalized response matrix describing the probability 
that a photon of energy E is detected in bin h, and A{E) is 
the detector effective area at energy E. Contrary to unfolded 
spectra produced by XSPEC (which are given by the model, 
rebinned to the output energy bins of the response matrix, 
multiplied by the data counts divided by the forward folded 
model counts, i.e., the ratio residuals), this definition produces 
a spectrum that is independent of the fitted model. In Fig. ^ 
we over plot the fitted model (using the internal resolution of 
the ‘ancillary response function’, i.e., the ‘arf’, used in the 
fitting process). The plotted residuals, however, are those ob¬ 
tained from a proper forward-folded fit. 

Although we previously have successfully fit the X- 
ray spectra with sophisticated Comptonization models 
(iNowak et al.l2003) . we obtain surprisingly good fits for nine 
of the ten radio/X-ray spectra using the following simple 
model (using the ISIS/XSPEC model definitions): absorp¬ 
tion (the phabs model, with Ah fixed to 6 x 10^* cm^) and a 
high energy, exponential cutoff (the highecut model) mul¬ 
tiplying a doubly broken power-law (the bkn2pow model, 
with the first break being in the far IR to optical regime, and 
the second break being constrained to the 9-12 keV regime) 
plus a gaussian line (with energy fixed at 6.4keV). Results 
are presented in Table 2. When considering just the X-ray 
spectra, a singly broken power-law fits all ten spectra, with 
better results than any of the Comptonization models that we 
have tried. The 12-200 keV power-law is typically seen to 
be harder than the 3-12 keV power-law by AF « -0.2. (In 
a future work, correlations of this spectral break with overall 
hardness will be presented for over 200 hard state spectra of 
Cyg X-1; Wilms et al., in prep.) 

The phenomenological power-law model was of course 
chosen because we are attempting to answer a phenomeno¬ 
logical question: do the extrapolated radio and X-ray spectra 
predict the amplitude of the IR flux, and the location of any IR 
break? The power-law provides the simplest model to extrap¬ 
olate. The additional power-law break in the 9-12keV band 
is required by the data. The energy band restriction in the fit¬ 
ting process was to avoid spurious local minima, caused by 
the break interfering with the gaussian line at low energy or 
with the highecut model at high energy. All but one ob¬ 
servation has break energy values, including 90% confidence 
level error bars, that fall well within this 9-12keV range. 

Phenomenologically speaking, the X-ray break is con- 
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Fig. 3.— Results of broken power-law fits to GX 339—4, showing the 
location of the break between the radio and soft X-ray power-law. Labels 
refer to Table 1, and are further described in the text. Here, and throughout 
the remaining figures, solid points will refer to obseiwations A (not shown in 
the above ) and B, which have ‘optically thin’ radio spectra. Dashed lines 
show the approximate integrated X-ray flux and IR spectral break energy 
previously observed in GX 339—4 (Corbel & Fender 2002). The dotted line 
is tx 


sistent with the expectatio ns from reflection models 
dMagdziarz & Zdziarskill9^ . However, as shown in Fig.|2] 
the identical broken power-law model described above pro¬ 
vides a very good phenomenological description even for 
Cyg X-1 hard state data with an extreme break. Although 
not mutually exclusive with the presence of reflection, such 
data require at least two broad-band continuum components 
in the X-ray. For example. Fig. |2l could be consistent with 
syn chrotron and s’ynchrotron self-Compton from jet mod¬ 
els dMarkoff & Nowakll2004 Markoff, Nowak, & Wilms, in 
prep.), or a strong disk plus Comptonization component from 
corona models. These issues will be explored in greater detail 
in a forthcoming work (Wilms et al., in prep.). 

2.2. Predicted Radio/X-ray Correlations 

Given only two radio points and an X-ray spectrum that 
is well-fit by a singly broken power-law, it is not surprising 
that the doubly broken power-law models work as well. But 
how do the fitted locations of the radio-to-soft X-ray break 
compare to the TR coincidence’ breaks, and how do the htted 
break locations scale with observed flux? In Fig. |3we show 
the fitted radio-to-X-ray break location as a function of 3- 
9 keV integrated flux. We also show in this figure the approx¬ 
imate integrated 3-9 keV flux a nd the IR break location f or the 
1997 observation discussed bv ICorbel & Fend^ (l2002h . For 
our GX 339-4 observations of comparable 3-9 keV flux, the 
doubly broken power-law models do indeed produce a break 
in the IR. Looking over the whole span of observed integrated 
X-ray fluxes, however, we see that the model fits presented 
here have predicted radio-to-X-ray breaks ranging all the way 
from the far IR to the blue end of the optical (and into the X- 
ray, if one also considers observation B [40108-01-04], which 
has an ‘optically thin’ radio spectrum). 

To assess the correlation of radio-to-X-ray break energy 
with integrated X-ray flux, we exclude observation B (40108- 
01-04), which has an optically thin radio spectrum, and per¬ 
form a regression analysis on the remaining eight data points. 
We weight the data uniformly, which is equivalent to assum- 


Fig. 4.— The flux of the radio power-law, integrated from zero energy to 
the spectral break between the radio and soft X-ray, vs. the 3-9 keV (PCA) 
flux. Labels refer to Table 1 (Observation A not shown), and are further 
described in the text. The dotted line is oc 


ing that intrinsic variations in any correlation dominate over 
statistical errors. (In the results discussed below, the derived 
regression slopes and errors encompass the values obtained 
if one weights the data by their error bars.) We find that the 
radio-to-X-ray break energy, in eV, scales with the 3-9 keV 
integrated flux as **. Does this agree with mod¬ 

els wherein the observed soft X-ray spectrum is the optically 
thin synchrotron emission from a jet? 

Using a scale invariance An s atz to descri be the jet 
physics (iHeinz & Sunvaevl l2Q03t [Heinz] 12004 see also 
iFalcke & Biermmnill995l and iMarkoff et alJl2003^ ■ we can 
predict the scaling between the integrated X-ray synchrotron 
flux and the radio-to-X-ray break frequency where the jet be¬ 
comes optically thin to synchrotron self-absorption. For sim¬ 
plicity, we make the assumption that, at the base of the jet, the 
fraction of particle pressure to magnetic pressure is indepen¬ 
dent of accretion rate/jet power. 

Given the above assumptions, the optically thick radio flux 
follows (iHeinz & Sunvae^2003^ : 


[ 2 p -{ p + 6 ) ocr + l 3 ] 
Tr OC B [f+4] 


( 2 ) 


where is the radio spectral index, p is the power-law index 
of the electron spectrum (at energies below any synchrotron 
cooling regi me), and B is the magnetic held at the base of 
the jet. From iHein using the same assumptions, the 

X-ray synchrotron hux integrated over a hxed energy band 
follows 

Fx cx , (3) 


where ofx is the X-ray spectral index. (We again are using the 
convention that a= 1 -F.) Combining these expresions, we 
hnd 

[2/?-(/?+6)Qr+13] 

Tr CX Fx [p+4][l+2p+3ax] . (4) 


If the X-ray band is unaffected by synchrotron cooling, the 
relation follows Tr oc values of 2.6 < p < 2, 

-0.8 < ttx = (1 -p)l2 < -0.5, and 0 < Ofr < 0.2. This scaling 
is consistent with previous observations, and is also consis- 
tent with prior descriptions of the jet mo del’s predictions (i.e., 
iCorhel et alJ^tMlMarkoff et alJ2003l) . 
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The lower break frequency, from an optically thick to 
optica lly thin spectrum is proportional to ilHeinz & ■Siinvaevl 
12001 

[6+p] 

i/bOcfiP^ . (5) 

We thus can write the dependence of z/b upon Fx as 

|6+Pl 2[p+6] 

Vh OCF^ ll+2p+3ax][4+p] oc [/’+4][p+5] ^ (g) 

where for the latter relation we have used ax = (1 - p)/2 from 
standard synchrotron theory. Again taking -0. 8 < ax < -0.5, 
we obtain i>\, oc to z/b c)c This prediction is flat¬ 

ter than the observed dependence of extrapolated break fre¬ 
quency upon X-ray flux. 

Eq. 0 represents the scaling of the actual location of the 
break from an optically thick to optically thin spectrum (as 
appropriate for the TR coincidence’ yielding the appropriate 
location of an observed break). Our observations instead yield 
the scaling of the inferred location of the break, which has 
an additional dependence upon the scalings of the spectral 
slopes, ar and a^. Specifically, if one assumes that there is 
an underlying radio/X-ray correlation given by oc F-J^, yet 
allows evolution of the radio and X-ray spectral slopes, one 
can show that for the inferred break, z/b, 


log 


( 

V^'o 


[(ar-ax)-l-Aar-Aax] ' x 


(1-/3) log 


-I-log 


^ Aar - log f 


\i^o 


V^'o 


Aax 


,(7) 


where z/q is the break energy measured at an X-ray flux of 
Fq, z/f, and z/x are the frequencies at which the radio and X- 
ray fluxes, respectively, are measured, and Aar and Aax are 
the changes in radio and X-ray spectral indices as the X-ray 
flux changes to Ex. (Note that small correction factors for the 
dependence of integrated X-ray flux upon ax, as opposed to 
solely power law normalization, have been omitted.) 

In the absence of evolution of the spectral indices, one 
expects the break frequency to scale with X-ray flux with 
a power of (1 -/3)/(ai-- ax ) ~ 0.35, as described above. 
Howe ver, as noted elsewhere (IWilms et al.ll999l:lNowak et al.l 
[in^, many of the observations of GX 339-4 show a soft¬ 
ening/hardening with increasing/decreasing X-ray flux. (For 
sources with optically thick radio spectra, notable exceptions 
to this trend are Obs IDs 40108-02-03 and 70109-01-02. 
A detailed discussion of flux/hardness trends can be found 
in iNowak et al.ll2002i ') Comparing Obs ID 40108-02-02 to 
20181-01-01, Attx « 0.06 for log(Ex/Eo) « 0.7. Combining 
this with ar-l-ax ~ 0.8 and log(z/x/z/o) ~ 4.9 means that the 
break energy should scale with with X-ray flux with a power 
of (l-/3-l-0.5)/0.8 « 1. This is in much closer agreement with 
the observed trends. 

Contrary to ax, there are no obvious trends for a depen¬ 
dence of tti upon X-ray flux. (This will be explored in a future 
work, with a larger set of observations from 2002 and 2004 
outbursts of GX 339-4; Corbel et al., in prep.) The mean 
value for all the optically thick radio data is (ai ) = 0.16, with 
a scatter of ±0.05. Combining this with log(z/r/z/o) « -3.5, 
and the fact that from the faintest to brightest observation is 
a factor of 10* we expect the variations in radio slope to 
account for correlation slope variations of up to ±0.1. 

We thus And that the observed correlation of inferred break 
energy is inconsistent with the expectations from the very 
simplest X-ray synchrotronjet model (i.e., the jet spectra scal¬ 
ing with input power, and no evolution of the radio or X-ray 
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Fig. 5.— 3-9keV (triangles) and 20-200keV (circles) X-ray flux (units 
of 10“* ergs cm s“*) from pointed RXTE observations of GX 339-4 vs. 8.6 
GHz radio flux (mjy). Labels refer to Table 1, and are further described in 
the text. Lines show power-laws of the form fx*, with b = 0.72 (dotted line) 
and b = 0.75 (dashed line). 
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slopes). The observed trends can be recovered, however, if 
the differences from the simplest expectations are primarily 
driven by evolution of the X-ray spectral slope, ax, with X- 
ray flux, with possible lesser contributions from variations of 
the radio spectral slope, a,. We shall discuss this point further 
in Cl 

A further interesting correlation is seen when one looks 
at the dependence of radio flux, integrated between zero 
energy and the break frequency, upon the integrated X-ray 
flux, as we show in Fig. Considering both in units of 
10“^ erg cm“^ s“*, and again weighting the data uniformly, we 
And Er = 0.04Ex* This implies that, even ignoring the 

energy of the particles and fields in the jet, the ‘jet radiation’ 
up to the inferred break is 3-10% of the integrated 3-9 keV 
X-ray flux, and is 0.5-5% of the 3-200 keV flux. 

Again, we can compare this result to expectations from the 
simplest jet models. Utilizing our previous estimates, we And 
for the product of TiVh'. 

[19-i-3p-ar(/?+6)] 

±;z/b oc //^L;>+4JLl+2p+3axJ (g) 

Again using the range of-0.8 < ax < -0.5 and 0 < ar < -0.2, 
this product should follow ^jz^, oc Ex*^* *®. As before, for the 
inferred value we might expect to increase the scaling expo¬ 
nent by 0.5 ±0.1, if we allow for an evolution of ar and ax 
with X-ray flux. 

3. RADIO/X-RAY CORRELATIONS IN GX 339-4 

One of the interesting aspects of the radio/X-ray flux corre¬ 
lation in GX 339-4 as identified bv iCorbel et al.N2003ft is that 
it appeared to be consistent^ in both slope and amplitude be¬ 
tween the 1997 and 1999 low/hard states, despite the interven¬ 
ing radio quiet, soft state outburst (see Fig.|3- The high lumi¬ 
nosity observations C and D [70109-01-02 and 40031-03-01], 
however, were not part of the original correlation presented 

^ There are two radio points, however, from the 1999 decline into quies¬ 
cence that deviate from the coirelation (Fig. 1^. Given_the_J^edegree of 
variability on week time scales in the quiescent state JKong_etalJ200^, and 
the fact that the magnitude of any time delay between the radio and X-ray 
flux variations is unknown, it was not clear how significant these deviations 
should be considered. 
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ASM Count Rate (cps) ASM Count Rate (cps) 

Fig. 6.— Left: 15 GHz Ryle radio flux (mJy) vs. Cyg X-1 daily mean ASM count rate (see text). Grey boxes highlight data associated with failed state 
transitions, and the beginning of a transition to a prolonged soft state. Right: 20-200 keV flux (units of 10“^ ergs cm s“^) from pointed RXTE observations, vs. 
Cyg X-1 daily mean ASM count rate (see text). 


hv ICorhel et^(l20nl . These two observations occurred in 
2002, 14 days apart, in a bright hard state as GX 339-4 was 
rising into an even brighter soft state. 

In Fig. 13 we plot the 8.6 GHz radio flux density"* vs. th e 
3-9 and 20-200keV X-ray fluxes. As in ICorhel et alH2003h . 
we have performed a regression fit for all the 1997 and 1999 
observations for which we could determine both a 3-9 keV 
{PCA) and 20-200keV (HEXTE) flux. (However, we also 
plot the 1999 low flux observations from ICorhel et iOlMl 
that have only PCA data.) For the soft X-rays we find 
jr = 9 9 ^^o. 73 ±o. 02^ X-rays we find = 

The slight differe nces between these correla- 
tion slopes and those presented bv ICorhel et al .1 il2003l) are due 
to our use of better calibrated RXTE response matrices (e spe- 
cially for the PCA). As discussed bv ICorhel et al.l(l2003h . the 
correlations apply equally well to the 1997 and the 1999 data, 
including the low flux 1999 data that were not formally part of 
the regression fit. The slopes between the soft an d hard X-rays 
are al so reasonably consistent with one another dCorhel et al.l 
12001 . 

The 2002 data, although only consisting of two points, de¬ 
viate from these previously viewed trends. For these data, the 
radio flux vs. the 3-9 keV flux correlation appears to have the 
expected slope, but more than a factor of two times lower am¬ 
plitude. The radio/hard X-ray correlation, however, appears 
to be less reduced in amplitude, but has a steeper slope, with 
J-r oc Tx* There is one major uncertainty in these results: 
the radio data associated with observations C and D were not 
strictly simultaneous with the X-ray data. We note, however, 
that the radio data, which falls below the previously observed 
trends, was taken after the X-ray data while the source was 
rising in flux on several day time scales. 

It is also possible that observations C and D, rather than 
being on a lower radio amplitude parallel track, are instead 
indicating a turnover in the correlation occuring at high X-ray 
flux. Such a turn over has been ob served for the correlation 
seen in Cyg X-1 (iGallo et al.ll2fl03l and Fig.|6j- If this latter 
hypothesis is correct, then given the observed 14.2 mJy flux 

Note that for observation C (70109-01-02), we have extrapolated the 
2.4 GHz radio flux to 8.64 GHz using the fitted power-law. 


for the highest flux observation, D, the deficit is a factor of 
2.3 for the 3-9 keV coiTelation (which predicts 32.4 mJy), but 
only a factor of 1.6 for the 20-200keV flux correlation (which 
predicts 22.7 mJy). Extrapolating the observed 2.4 GHz radio 
flux for observation C to 8.64GHz yields 5.9 ±0.9mJy. This 
is a factor of 2 );q j below the value predicted by the radio/3- 
9keV flux correlation, and a factor of 2.2)^-3 below the value 
predicted by the radio/20-200 keV flux correlation. 

A further interesting fact to note also comes from the ra¬ 
dio/hard X-ray correlation. Observations A and B (40108-01- 
03 and 40108-01-04) were the first observations taken after 
the 1999 return to the radio-loud hard state, and both have op¬ 
tically thin (i.e., < 0) radio spectra. Although they have 

nearly identical 3-9 keV integrated fluxes, they have different 
radio fluxes and they straddle the radio/soft X-ray flux cor¬ 
relation shown in Fig. |3 However, these points actually lie 
along the radio/hard X-ray flux correlation, as they have sig¬ 
nificantly different 20-200keV fluxes. This might argue that 
the radio/hard X-ray flux is the more fundamental relation¬ 
ship. 

4. RADIO/X-RAY CORRELATIONS IN CYG X-1 

The radio/X-ray correlations in GX 339-4 present us with 
two, not necessarily mutually exclusive, possibilities. Either 
the radio/X-ray correlation is more fundamentally tied to the 
hard X-rays, and/or different instances of the hard state can 
present the correlation with different amplitudes. To further 
explore these pos sibilities, we turn to radio/X-ray observa - 
tions of Cyg X-1 l|Pottschmidt et al J|2()( )3HGallo et al. 11200^ . 
As described in iPottschmidt et alJ J2003h . the radio data are 
15 GHz observations performed at the Ryle Telescope, Cam¬ 
bridge (UK). These are single channel observations, so a radio 
spectral slope cannot be determined. Most of these observa¬ 
tions have occurred simultaneously with pointed RXTE ob¬ 
servations llPottschmidt et al.ll2003l) . and nearly all have very 
good contemporaneous coverage by the RXTE All Sky Mon¬ 
itor (ASM). The ASM provides covera ge in the 1.5-12keV 
energy band JRemillard & Levine! 1997^ ■ and therefore allows 
us to assess the radio/soft X-ray correlation in Cyg X-1. 

In Eig.|3we plot the daily average ASM count rate vs. the 
daily average 15 GHz flux. Here we only include days for 
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which there are at least 25 ASM data points (representing 70 
second d wells) with a reduced < 1.5 for the ASM solu¬ 
tion (see lRemillard & Levine^lT997^ . Ranging from approxi¬ 
mately 10-50cps in the ASM there is a clear log-linear cor¬ 
relation between the radio flux and the ASM count rate. As 
for GX 339-4, the radio flux rises more slowly than the ASM 
count rate (Ti- scales approximately as the 0.8 power of the 
ASM count rate). Cyg X-1, however, shows much more scat¬ 
ter in the amplitude of the correlation than does GX 339-4. 
This scatter is not dominated by va riations related to the o r- 
bital phase of the binary system (cT. lBrocksonn et al.ll999h . 

As noted bv iGallo et alJ(l2003h . there is a sharp roll-over for 
higher ASM count rates. However, one can clearly discern on 
the shoulder of this roll-over (i.e., the upper right corner of 
Fig-EJ four ‘spokes’, consisting of 2-5 data points each. In 
these spokes, the radio/X-ray correlation appears to hold to 
high count rates. We have confirmed^ that each of these times 
are associ ated with ‘failed transitio ns’ to the soft state, as de¬ 
scribed bv iPottschmidt et al.l il2003h . except for the lowest am¬ 
plitude of these spokes, which occurs immediately preceding 
a prolonged soft state outburst. In many ways these failed 
state transitions and the early stage of a soft state outburst for 
Cyg X-1 are reminiscent of the data presented here for the 
2002 GX 339-4 hard state. As for those latter GX 339-4 
data, the Cyg X-1 spectra remain hard, although not as hard 
as low-luminosity hard states, to very high flux levels. 

In Fig. 1^ we also plot the daily average ASM count rate 
vs. the 20-200keV flux from our pointe d RXTE observa¬ 
tions taken during the same 24 hour period jPottschmidt et alJ 
l2003h . Here, however, we only require 3 ASM data points per 
average, for a resulting 85 observations. We see that hard X- 
ray/ASM correlation traces a similar pattern to the radio/ASM 
correlation. Indeed, when we plot the hard X-ray flux vs. 
the daily average radio flux (representing 129 separate hard 
X-ray pointed observations), we obtain a log-linear re lation- 
ship, as shown in Fig.0(see also lGleissner et al.l2flfl4l) . (Two 
high radio flux, but low X-ray flux, data points are seen in 
the upper left of this figure; these represent intra-day hard X- 
ray variations.) Returning to the question that we posed in 
of whether we were seeing evidence for parallel tracks in 
the radio/X-ray correlation, or whether we were seeing the 
correlation being more fundamental to the hard X-ray band, 
in Cyg X-1 the answer seems to be that both effects are dis¬ 
cernible, but fundamentally the radio flux density does appear 
to be tied to the hard X-ray emission. 

5. SUMMARY 

In this work we have considered ten simultaneous or near 
simultaneous RXTE/radio observations of GX 339-4, and 
over one hundred RXTE/radio observations of Cyg X-1. We 
have fit the former spectra with a very simple, but remark¬ 
ably successful, phenomenological model consisting of a dou¬ 
bly broken power-law with an exponential roll-over plus a 
gaussian line. For GX 339-4, the break between the radio 
and soft X-ray power-laws occurs in the IR to optical ra nge, 
in agreement with the p rior work of iMotch et alJ (119851) and 
ICorbel & Fended (l2002^ (i.e., the ‘IR coincidence’). In con¬ 
trast to these prior works, we have fit the X-ray data in ‘de¬ 
tector space’ and provided a quantitative assessment of the 

^ We have used the ‘vwhere’ program INoblel2()05ll which allows one to 
interactively filter one set of data, e.g., a color-color or color-intensity dia¬ 
gram, and then apply those filters to the same observations visualized in a 
different way, e.g., a time-intensity diagram. (The program is available at 
http://space.mit.edu/cxc/software/slang/modules /1 



20-200 keV Flux (10 ^ ergs cm ® s ‘) 

Fig. 7.— 20-200keV flux (units of 10“^ ergs cm s“*) vs. tlie daily average 
15 GHz Ryle radio flux (mJy) for pointed observations of Cyg X-1. 


extrapolated break location. 

Is the ‘IR coincidence’ just a coincidence? The spectral 
fits to the GX 339-4 data suggest that the answer is ‘possi¬ 
bly not’. All of the data with optically thick radio spectra 
are extremely well-fit by doubly broken power-laws, with the 
break being in the IR regime. Although the scaling of this 
break frequency with X-ray flux does not agree with the pre¬ 
dictions of the simplest X-ray synchrotron jet models, if one 
allows for softening of the X-ray spectrum with increasing 
X-ray flux, then the jet model predictions agree with the ob¬ 
served correlations. There are several possibilities why such 
a softening could occur. In terms of the jet model, for exam¬ 
ple, there could be cooling of the electron spectrum leading 
to a steepening in the X-rays. A second possibility is that the 
soft X-rays are being contaminated by excess emission from 
a disk component that becomes more prominent with increas¬ 
ing flux. What is clear is that in order for jet models to explain 
both the scaling of oc E, and i/j, oc , they must be 
more sophisticated than the very simplest considerations. Fits 
to some of these data sets with such models, which include 
disk emission, synchrotron self-Compton emission, etc., are 
currently being considered in detail, and will be discussed in 
a future work (Markoff, Nowak, & Wilms, in prep.). The re¬ 
sults presented here, however, are consistent with the notion 
that at least some fraction of the observed soft X-rays may 
be attributable to emission from the jet, as opposed to disk or 
corona. 

On the other hand, other facts suggest that the answer to the 
question of the ‘IR coincidence’ being just that is ‘maybe’. 
We have tentative evidence in GX 339-4, and firmer evi¬ 
dence in the Cyg X-1 failed state transitions and soft state 
transition, that the correlation between radio flux and inte¬ 
grated X-ray flux can take on different amplitudes during dif¬ 
ferent hard state episodes. There is also evidence in Cyg X-1 
that the radio/X-ray correlation is more fundamental to the 
hard X-ray band. In jet models, this band, which essen¬ 
tially encompasses the third, highest energy, power-law com¬ 
ponent in our model fits (and also encompasses the exponen¬ 
tial cutoff), is possibly attributable to the synchrotron self- 
Compton (SSC) emission from the base of th e jet (see, e.g., 
iMarkoff et al.ll2003UMarkoff & Nowakll2004|) . It is therefore 
quite reasonable to expect a strong coupling between the radio 
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and hard X-ray flux; however, these models are more com¬ 
plex than simple pure synchrotron mode ls, and are only now 
beginning to be explored quantitatively dMarkoff et al.l(20()^ 
iMarkoff & Nowakll2004l) . It is also worth noting that such 
radio/X-ray couplings can be expected in Comp ton corona 
models (e.g., lMeieTl200lHMerloni & Fabianl200^ . 

More problematic for pure jet models are GX 339-4 obser¬ 
vations A and B (40108-01-03,40108-01-04), both of which 
occurred very shortly after the 1999 return to the hard state. 
They have ‘typical’ hard state properti es in terms of spec¬ 
tra and variability JNowak et al.l 1200^ . but have optically 
thin radio spectra that do not neatly extrapolate into the ob¬ 
served soft X-ray spectra. They satisfy the radio flux den¬ 
sity/integrated X-ray flux correlation; however, they agree 
better with the radio/hard X-ray flux correlation, for which 
these two observations are cleanly differentiated from one an¬ 
other (Fig.|3l. It is possible that they represent an early, tran¬ 
sient phase of the jet as it forms. It then remains a theoreti¬ 
cal question for the jet model whether it can describe a sce¬ 
nario where the basic radio/X-ray flux correlation and ‘typi¬ 
cal’ steady-state X-ray spectra hold, while the radio has not 
yet settled into a steady ‘optically thick’ state. 

The results presented here suggest, at the very least, some 
obvious observational strategies to arrive at a more definitive 
answer to the question that we have posed. Given the break 
energy correlations, it would be extremely useful to have not 
only a radio amplitude for each X-ray observation, but also a 
radio slope. Furthermore, the inferred break for the brightest 
observation, D (40031-03-01), occurs in the blue end of the 
optical. Thus, ideally multi-wavelength observations would 
consist of radio, broad band X-ray, and IR through optical 
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cover age (see, for example, iHvnes et alJl2003t iMalzac et alJ 
l2(M which suggest the possibility of jet synchrotron radi¬ 
ation extending into the optical regime for the hard state of 
XTE J11183-48). This is an admittedly difficult task, but BHC 
are demonstrating via spectral correlations that all these en¬ 
ergy regimes are fundamentally related to activity near the 
central engine, i.e., the black hole, of the system. 

Finally, it is important to obtain multi-wavelength observa¬ 
tions of multiple episodes of each of the spectral states. For 
example, if there are indeed ‘parallel tracks’ in the radio/X- 
ray correlations, it would be interesting to determine whether 
the amplitude of the radio/X-ray correlation is related to the 
flux at which the outbursting source transits from the low/hard 
to high/soft state. We note that the low amplitude radio/X- 
ray correlation of the 2002 GX 339-4 outburst was associated 
with a very high flux level for the hard-to-soft state transition. 
The low amplitude ‘spokes’ in the Cyg X-1 radio/X-ray cor¬ 
relation might be similar in that they are associated with failed 
state transitions. If such observations can be made with more 
quantitative detail, we will have vital clues to determining the 
relative contributions of coronae and jets, and the coupling be¬ 
tween these two components, for black hole binary systems. 
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TABLE 1 

GX 339-4 Observed Fluxes (la- confidence). 


ObsID 


y.m.d 

0.84 GHz 

1.38 GHz 

2.40 GHz 

4.80 GHz 

8.64 GHz 

3-9 keV 

9-20 keV 

20-100 keV 

100-200 keV 






(mJy) 




( 10 -'* 

ergs cm~^ s“^) 


20181-01-01 


1997.02.03 

7 0^-^ 




^•^- 0.1 

0.90 

0.92 

2.98 

1.27 

20181-01-02 


1997.02.10 

'^•'^-0.7 





0.80 

0.82 

2.67 

1.11 

20181-01-03 


1997.02.17 

5 6^-'^ 


o.u_o 2 



0.77 

0.78 

2.63 

1.06 

40108-01-03 

A 

1999.02.12 




^ 04+0.08 
^•■^^-0.08 

^•‘^^-0.08 

0.41 

0.39 

1.18 

0.48 

40108-01-04 

B 

1999.03.03 




5 Q7+0.06 
^•^^-0.06 

c 74+0.06 
' -0.06 

0.42 

0.42 

1.51 

0.66 

40108-02-01 


1999.04.02 




A -7«:+0.06 
^•'•^-0.06 

5 1Q-^0.06 
^•^^-0.06 

0.42 

0.44 

1.60 

0.97 

40108-02-02 


1999.04.22 




0 QO+0.06 
^•^^-0.06 

■^•^^-0.06 

0.19 

0.21 

0.74 

0.43 

40108-02-03 


1999.05.14 




1 oc+U.Ub 
^•"^'^-0.06 

^ -^^-0.06 

0.07 

0.07 

0.21 

0.14 

70109-01-02 

C 

2002.04.03 


A 071+0.20 
^•°-^-0.20 

c 10+0.11 
J.iO_0 11 



1.30 

1.47 

5.34 

2.09 

40031-03-01 

D 

2002.04.18 




-0.06 

i^t.iO-Oo? 

5.10 

4.97 

12.3 

2.42 


Note. — X-ray flux errors are taken to be 3.5%, the average of the fitted normalization difference between PCA and HEXTE, since this exceeds the statistical uncertainties. 


TABLE2 

GX 339-4 Spectral Fits (90% Confidence Level errors). 


Obs ID 


Ecut 

Efoid 

Qr 

Eb-r 

E soft 

Eb-x 

Ehard 

4^1ine 

•^line 

Con. 

X^/DoF 



(keV) 

(keV) 


(eV) 


(keV) 


(10-3 ) 

(keV) 



20181-01-01 


37+7 

^'-6 

15512 

0 11+0-08 
^•^^-0.07 

0.801:22 

1 69+0.01 
^•^^- 0.01 

11 3+0-6 

1 49+0.03 
^•^^-0.02 

2 1+0-2 

^••^-0.3 

0 6+0-1 
^•^-0.1 

0 97+0.01 
'^•^'-0.02 

110/108 

20181-01-02 


42+5 

1401® 

Q 17+0.10 

^-0.09 

0.37l:“ 

1 6;s+ooi 
^•°°-o.oo 

1 1 n+0.6 

l^-^-0.5 

1 40+0.02 
^•^°-0.02 

-1 O+0.3 

l-°-0.2 

0 7+0-1 
'-0.2 

0.951:111 

162/114 

20181-01-03 


53!r 

133?1 

^•^‘^-0.04 

1-391:1 

1 618+0.01 
^•°°-o.oo 

10 7+0-6 

1 ci-bO.02 
^•'^^-0.02 

1 7+0-2 

1- '-0.2 

0 6+0-1 
^•‘^-0.1 

0.98l°2 

90/90 

40108-01-03 

A 

441, 

1651? 



1 8n+ooi 
E^^-o.oi 

10 2+1-0 

1 59+0.02 
^•-^^-0.03 

-1 2+0.1 

l-^-0.2 

0 7+0-1 
^•'-0.1 

1 03+0-01 
+ -D.7-0 01 

64/77 

40108-01-04 

B 

69!1 

1381* 

^•^^-0.05 

13-51*0 

1 69+0.01 

^•^^-0.01 

10 

1U.O_0.7 

1 C9+O.OI 
^•'-’■^-0.02 

^•^-0.1 

u.O -02 

1 03^-02 
^•^■^-0.02 

106/103 

40108-02-01 


321‘ 

34511 

Q 19+0.06 
^•^^-0.05 

0-441:1 

1 67+001 

10.51:® 

1 45+0.03 
^•^'^-0.02 

0 7+0-2 
^•^-0.1 

0 5+0-2 
^•^-0.2 

0 97+0.02 
^•^'-0.02 

115/93 

40108-02-02 


44+33 

^^-18 

304122 

^•^^-0.08 

0 12+0-18 
^•^^-0.06 

^•^•^-0.01 

10.71:2 

1 47+0.03 
^•^'-0.05 

^•^-0.1 

u.j-o_2 

0 95+0.03 
^•^‘'-0.03 

83/67 

40108-02-03 




0 24+^' 

0 07+0-16 
^•^'-0.04 

1 6;s+ooi 
ED6 _oo2 

9 0+2- 
^•^-0.0 

1 62+0-02 
^•^^-0.05 

o-i1:i 

0 6+0-2 
-0.3 

0 95+0.07 
^•^'^-0.08 

26/44 

70109-01-02 

C 

34l 

II 2 I 

0 1 1+0-14 
^•1^-0.11 

Q 95+14.3 
^•^^-0.81 

1 ^1+0.01 

10 3+0-4 

1'0-.7_o4 

1 97+0.01 
^•■^'-0.02 

2 - 21 ? 

0-41:1 

0 99+0.01 
^•^^-0.01 

180/148 

40031-03-01 

D 

26l 

69l 

0 0 

d 

9 OO+0.69 
-0.41 

1 79+0.01 
^•'•^-0.01 

11 l'^-^ 

-1 C9+O.O2 

^•--^^-0.03 

17.71:° 

0.811 

0 99+0.02 
^•^^-0.01 

143/130 


Note. — 20181 data are from the 1997 hard state (Wilms et al. 1999; Nowak, Wilms, & Dove 2002), 40108 data are from the 1999 hard state (post-1998 soft 
state; Nowak, Wilms, & Dove 2002), 70109 and 40031 data are from the 2002 very luminous hard state (post a 3 year duration quiescent state; Homan et al. 2004). 
Parameters in italics were held fixed at that value for the fits. Note that the radio spectral index, ar, is defined to be positive for ‘optically thick’ spectra, and is related 
to the usual X-ray astronomy definition of the photon index by ar = 1 -Lr. 









